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In this lecture, we will focus on examples of electromagnetic wave propagation in nonlinear optical 
media, by applying the forms of Maxwell's equations as obtained in the eighth lecture to a set of 
particular nonlinear interactions as described by the previously formulated nonlinear susceptibility 
formalism. 

The outline for this lecture is: 

• General process for solving problems in nonlinear optics 

• Second harmonic generation (SHG) 

• Optical Kerr-effect 

General process for solving problems in nonlinear optics 

The typical steps in the process of solving a theoretical problem in nonlinear optics typically 
involve: 

J define the optical interaction of interest 
1^ (identifying the susceptibility) 

define in which medium the interaction take place 1 
(identify crystallographic point symmetry group) j 

j consider eventual additional symmetries and constraints 1 
1^ (e. g. intrinsic, overall, or Kleinman symmetries) j 

J construct the polarization density 1 
1^ (the Butcher and Cotter convention) j 

formulate the proper wave equation 
(e. g. taking dispersion or diffraction into account) 

J formulate the proper boundary 
1^ conditions for the wave equation 

J solve the wave equation under 1 
1^ the boundary conditions j 
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Formulation of the exercises in this lecture 

In order to illustrate the scheme as previously outlined, the following exercises serve as to give 
the connection between the susceptibilities, as extensively analysed from a quantum-mechanical 
basis in earlier lectures of this course, and the wave equation, derived from Maxwell's equations of 
motion for electromagnetic fields. 

Exercise 1. {Second harmonic generation in negative uniaxial media) Consider a continuous pump 
wave at angular frequency w, initially polarized in the y-direction and propagating in the positive 
a;-direction of a negative uniaxial crystal of crystallographic point symmetry group 3m. (Examples 
of crystals belonging to this class: beta-BaB204/BBO, LiNbOs.) 

la. Formulate the polarization density of the medium for the pump and second harmonic wave. 

lb. Formulate the system of equations of motion for the electromagnetic fields. 

Ic. Assuming no second harmonic signal present at the input, solve the equations of motion for 
the second harmonic field, using the non-depleted pump approximation, and derive an expression 
for the conversion efficiency of the second harmonic generation. 

Exercise 2. [Optical Kcrr-effect - continuous wave case) In this setup, a monochromatic optical 
wave is propagating in the positive ^-direction of an isotropic optical Kerr-medium. 

2a. Formulate the polarization density of the medium for a wave polarized in the a;y-plane. 

2b. Formulate the polarization density of the medium for a wave polarized in the x-dircction. 

2c. Formulate the wave equation for continuous wave propagation in optical Kerr-media. The 
continuous wave is x-polarized and propagates in the positive ^-direction. 

2d. For lossless media, solve the wave equation and give an expression for the nonlinear, 
intensity-dependent refractive index n = no + n2\E^\^ . 
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Second harmonic generation 

The optical interaction 

In the case of second harmonic generation (SHG) , two photons at angular frequency w combine to 
a photon at twice the angular frequency, 

hv + hw ^ h(2oj). 

This interaction is for the second harmonic wave (at angular frequency a;)described by the second 
order susceptibility 

x|fi/3(-Wa;w,a;), 

where = 2a; is the generated second harmonic frequency of the light. 

Symmetries of the medium 

In this example we consider second harmonic generation in trigonal media of crystallographic point 
symmetry group 3m. (Example: LiNbOa) 



Ei(r,t) = Re[E^ exp(-iwt)] 

= eyRe[Ey exp{—iujt)] 



ET(r, t) = Re[E^ exp{-iojt)] + Re[E2u; exp{-i2u>t)] 




Crystal frame (a;, y, z) 



point symmetry group 3m 
(e. g. LiNbOa) 



Figure 1. The setup for optical second harmonic generation in LiNbOs. 

For this point symmetry group, the nonzero tensor elements of the first order susceptibility are 
(for example according to Table A3.1 in The Elements of Nonlinear Optics) 



(1) 



A,ZZ ' 



which gives the ordinary refractive indices 

n^{u) = Uyiu) = [1 +xS(-'^;'^)]^^^ = »^o(w) 
for waves components polarized in the x- or y-dircctions, and the extraordinary refractive index 

for the wave component polarized in the z-direction. Since we here are considering a negatively 
uniaxial crystal (see Butcher and Cotter, p. 214), these refractive indices satisfy the inequality 

n^ioj) < no(w). 

The nonzero tensor elements of the second order susceptibility are (for example according to 
Table A3. 2 in The Elements of Nonlinear Optics) 



y(2) =y(2) =y(2) 
Ayxx Axyx A,xxy 



.y(2) 

Ayyyi 



y(2) =T^(2) 



y(2) =^(2) 



y(2) 

y(2) =^(2) 



(1) 
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Additional symmetries 

Intrinsic permutation symmetry for the case of second harmonic generation gives 

which reduces the second order susceptibility in Eq. (1) to a set of 11 tensor elements, of which only 
4 are independent. (We recall that the intrinsic permutation symmetry is always applicable, as 
being a consequence of the symmetrization described in lectures two and five.) Whenever Kleinman 
symmetry holds, the susceptibility is in addition symmetric under any permutation of the indices, 
which hence gives the additional relation 

i. e. reducing the second order susceptibility to a set of 11 tensor elements, of which only 3 are 
independent. 

To summarize, the set of nonzero tensor elements describing second harmonic generation under 
Kleinman symmetry is 



(2) =^(2) ^^(2) =_^(2) ^(2) 

(2) 



Ayxcc Axj/cc Axxy '^yyy'' ■^zzz'i 



y(2) =yC^) =yC^) =y(2) =y(2) =y(2) 

A.ZXX ^zyy ^yyz A.xxz ^yzy A.xzx' 

For the pump field at angular frequency w, the relevant susceptibility describing the interaction 
with the second harmonic wave is^ 

x2!;3(-w;2w,-a;). 

For an arbitrary frequency argument, this is the proper form of the susceptibility to use for the 
fundamental field, and this form generally differ from that of the susceptibilities for the second har- 
monic field. However, whenever Kleinman symmetry holds, the susceptibility for the fundamental 
field can be cast into the same parameters as for the second harmonic field, since 

X^a^(— w; 2a;, —a;) = {Apply overall permutation symmetry} 

= xS/3(2w; -a;, -a;) 

= {Apply Kleinman symmetry} 

= x2/3(2a;; -a;, -a;) 

= {Apply reality condition [B.&C.Eq. (2.43)]} 

= [x2/3(-2^;'^''^r 

= x2/3(-2w;a;,a;). 

Hence the second order interaction is described by the same set of tensor elements for the funda- 
mental as well as the second harmonic optical wave whenever Kleinman symmetry applies. 



Keep in mind that in the convention of Butcher and Cotter, the frequency arguments to the 
right of the semicolon may be writen in arbitrary order, hence we may in an equal description 
instead use 

xEU-'^; -'^,2cj) 

for the description of the second order interaction between light and matter. 
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The polarization density 

Following the convention of Butcher and Cotter,^ the degeneracy factor for the second harmonic 
signal at 2uj is 

ii:(-2a;;a;,a;) = 2'+™->, 



where 



p = {the number of distinct permutations of = 1, 

n = {the order of the nonlinearity} ~ 2, 

m = {the number of angular frequencies LOh that are zero} = 0, 
1, if 2co ^ 0, 



0, otherwise. 



1. e. 



1, 



K{-2lj;lo,uj) = 2^+°"^ x 1 = 1/2 



For the fundamental optical field at u), one might be mislead to assume that since the second 

order interaction for this field is described by an identical set of tensor elements as for the second 
harmonic wave, the degeneracy factor must also be identical to the previously derived one. This 
is, however, a very wrong assumption, and one can easily verify that the proper degeneracy factor 
for the fundamental field instead is given as 

ii:(-w;2w,-w) = 2'+'"->, 



where 



p = {the number of distinct permutations of 2co, —lo} = 2, 
n = {the order of the nonlinearity} = 2, 

m = {the number of angular frequencies Uk that are zero} = 0, 

1, if uj^ 0, 



0, otherwise. 



= 1, 



1. e. 



K{-uj; 2uj, -w) = 2^+°-^ X 2 = 1. 
The general second harmonic polarization density of the medium is hence given as 

[P^i = [P^'j]. = £0 K{-2u,;u>,u;)x%{-2u;;u,,w) E^E^ 
' ^ ' 

= 5xfaV-2a';a-,c.) 

= {eo/2)[x^^lEZE- + xflElEl + x^^E^E^J 

= {^MxfL{EZEZ + ElEl) + xflE^EZ], 
Pr^]. = {eo/2UXKEZ + X^^ElEl + xf^ElEZ + xflyE^El] 

= {eo/2){x^^l{EZEZ - ElEl) + xWME^ + E^El)], 
[P^^]. = {eo/2)[x^XKEl + x'ilElEZ + x^^Ie^E^ + X^^Le^EZ] 

= {eo/2)[xfl{EZEl + ElEZ) + xfl{EZEZ + E^EZ)], 



^ See course material on the Butcher and Cotter convention handed out during the third lecture. 
Notice that for the first order polarization density, one at optical frequencies aiwajs has the trivial 
degeneracy factor 

K{-2iO-uj) = 2'+'"-> = 2i+°-i X 1 = 1. 
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while the general polarization density at the angular frequency of the pump field becomes^ 
[P£'''^']. = [P£'^]. = ^0 K{-uj- 2a;, -uj)x%{-uj; 2a;, -a;) E^^E^^ 

= eolx'^'liELEZ* + El^El*) + xf}AEl*l 
[P^^^]. = e,[xfLE^2u,EZ* + x'^JyELEZ* + xf^AE^J + xflyE^El*] 

= ^ob6yLiELEZ* - El^El*) + xfliEl^E^J + E^^Ey*)], 
i^n^ = ^o[x%E^^^El* + xflEl^EZ* + x^^lE^EZ* + X^^^AE^*] 

= ^oix^yliELEZ* + El^EZ*) + xfLiE^EZ* + E^.^EZ*)]- 

For a pump wave polarized in the yz-pla,ne of the crystal frame, the polarization density of the 
medium hence becomes 

[P^^]. = [eomix^^lEyEy+x'-^lE^E:,], 

[Pr\ = {^omi-X^X^^Ey + x'-^liE^EZ + EZEl)], 

[Pr]x=0, 

and 

rp(NL)i -£^[^(2) Ey Ey* + Y^^^ E^ E^*\ 

[Pif = eoi-X^^J AEy* + xfLiElA + E^^^El*)], 

[pir"^].=o. 

The wave equation 

Strictly speaking, the previously formulated polarization density gives a coupled system between 

the polarization states of both the fundamental and second harmonic waves, since both the y- 
and z-components of the polarization densities at u) and 2a; contain components of all other field 
components. However, for simplicity we will here restrict the continued analysis to the case of a 

(2) (2) 

y-polarized input pump wave, which through the Xzyy = Xzxx elements give rise to a ^-polarized 
second harmonic frequency component at 2a;. 

The electric fields of the fundamental and second harmonic optical waves are for the forward 
propagating configuration expressed in their spatial envelopes A^^ and A2u> as 

E^(a;) = ByAKx) exp{ik^^x), k^j^ = wn^^Jc = uino{u))/c 
^2uj{x) = ezA^^{x) exp{ik2o,,x), k2oj^ = 2u;n2ojJc = 2ujn^{2uj)/c 

Using the above separation of the natural, spatial oscillation of the light, in the infinite plane wave 
approximation and by using the slowly varying envelope approximation, the wave equation for the 
envelope of the second harmonic optical field becomes (see Eq. (6) in the notes from lecture eight) 

= ^^(|~~7^ yX?ix^S^^exp(2iA;i^^a;)exp(-ifc2a,,a;) 

(2) 

= i ^y2 exp[i(2fca,^ - k2oj^)x], 

^ri2ujz ^ 



^ Keep in mind that a negative frequency argument to the right of the semicolon in the suscep- 
tibility is to be associated with the complex conjugate of the respective electric field; see Butcher 
and Cotter, section 2.3.2. 
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while for the fundamental wave, 



dx 2k^ 

,2 



2{ujn^Jc) 



(2) 

These equations can hence be summarized by the coupled system 



(2) 



l- 



Ayexp{iAkx), (3a) 



^=.^AL^rexp(-.AU). (3b) 

where 

Afc = 2k^ - 

= 2ujn^^/c - 2Lon2^Jc 
= {2uj/c-}{n^^ - ujn2u,J 

is the so-called phase mismatch between the pump and second harmonic wave. 
Boundary conditions 

Here the boundary conditions are simply that no second harmonic signal is present at the input, 

^L(o) = o, 

together with a known input field at the fundamental frequency, 

Al{0) = {known}. 

Solving the wave equation 

Considering a nonzero Afc, the conversion efficiency is regularly quite small, and one may approx- 
imately take the spatial distribution of the pump wave to be constant, A'^{x) ~ A^{0). Using this 
approximation*, and by applying the initial condition ^2w(0) = of the second harmonic signal, 
one finds 







L (2) 

2(Q) exp(iAA;a;) dx 







2n2uj^c 



2 



A^^(O) — [exp(iAfcL)-l] 



= {Use [exp{iAkL) - 1]/Ak = iLexp{iAkL/2) sinc(AfcL/2)} 

.(2) 

-Al^{0) exp{iAkL/2) sinc(AA:L/2). 



2n2t^,c 



* For an outline of the method of solving the coupled system (1) exactly in terms of Jacobian 
elliptic functions (thus allowing for a depleted pump as well), see J. A. Armstrong, N. Bloembergen, 
J. Ducuing, and P. S. Pershan, Phys. Rev. 127, 1918-1939, (1962). 
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0.25 - 



_s 

^ 0.2- 



0.1 - 



0.05 - 




Figure 2. Conversion efficiency l2u>{L) / lu;{0) as function of normafized crystal length AkL/2. The 

conversion efficiency is in the phase mismatched case (Afc 7^ 0) a periodic function, with period 
2Lc, with Lc = 7r/Afc being the coherence length. 

In terms if the light intensities of the waves, one after a propagation distance x = L hence has the 
second harmonic signal with intensity /2w [L) expressed in terms of the input intensity Z^) as 

1 



hu.{L) = -socn2^MLiL)\ 

.'-^Xzxx-Li 



= ■^eocn2u>. 



■ £0 



AyJ{Q) exp(iAA:L/2) sinc(AA:L/2) 



Use \Aim^ 



2n2aj,C 

lxiinA^(0)rsinc2(AfcL/2) 



= £0 



IX 



(2) ,2 



21 

£ocn, 



sinc2(Afci/2) 



2£oc3 



i. e. with the conversion efficiency 

/2.(i) Lo'L^ |xixl(-2a;;a;,a;)|2 



(0) sinc2(AA;L/2), 



7^(0) 2£oc3 



1^(0) sinc2(AfcL/2). 
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1.52 - 



1.5' 1 1 1 

0.5 1 1.5 

vacuum wavelength 

Figure 3. Ordinary and extraordinary refractive indices of a negative uniaxial crystal as function 
of vacuum wavelength of the light, in the case of normal dispersion. Phase matching between the 
pump and second harmonic wave is obtained whenever n^,^ = no{ui) = tie (2a;) = n^u^- 

0.4 1 1 1 \ n \ 1 — \ 1 — I 1 n \ 1 1 



0.35 - 



0.3 - 



0.25 - 

_3 



^ 0.2- 




AkL72 

Figure 4. Conversion efficiency /2aj(-^)/^w(0) as function of normalized crystal length AA;L/2 when 
the material properties are periodically reversed, with a half-period of Lc- 
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Ayyzz 


Azzyy 


- y(3) - 


y(3) 


- y(^) 

A.xxyy 


= y' 

Ai 


y(3) 
^yzyz 


= yf-'^' 

A.zyzy 


- y(3) - 


y(3) 


- y(^^ 

Axyccy 


= y' 

Ai 


y(3) 
A.yzzy 


-y(3) 
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- y(^^ - 

A.ZXXZ 


y(3) 
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= y' 
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Optical Kerr-efFect - Field corrected refractive index 

As a start, we assume a monochromatic optical wave (containing forward and/or backward prop- 
agating components) polarized in the a;y-plane, 

E(^,t) = Re[Ea,(2)exp(-iwt)] e R^ 

with all spatial variation of the field contained in 

The optical interaction 

Optical Kerr-effect is in isotropic media described by the third order susceptibility^ 

(3) / ^ 

Symmetries of the medium 

The general set of nonzero components of Xfia^-y ^^"^ isotropic media are from Appendix A3. 3 of 
Butcher and Cotters book given as 



(3) 

yyxx 

yyxyx 

.(3) 
y-yxxy 

with 

y(3) =y(3) ,^(3) , (3) 
A.xxxx A.xxyy ' A.xyxy ' A-xyyx^ 

i. e. a general set of 21 elements, of which only 3 are independent. 
Additional symmetries 

By applying the intrinsic permutation symmetry in the middle indices for optical Kerr-effect, one 
generally has 

which hence slightly reduce the set (4) to still 21 nonzero elements, but of which now only two are 
independent. For a beam polarized in the xy-plane, the elements of interest are only those which 
only contain a; or y in the indices, i. e. the subset 

y(3) ^ ^(3) (3) ^ (3) ^ f intr. perm. symm. 1 ^ (3) ^ (3) (3) ^ (3) 

A.XXXX Ayyyy^ Axxyy Ayyxx | (a, a;) ;=i j -^xyxy Ayxyx' Axyyx -^yxxy^ 

with 

(3) ^ (3) , (3) , (3) 
■\xxxx Axxyy ' Axyxy ' ^xyyx^ 

i. e. a set of eight elements, of which only two are independent. 

^ Again, keep in mind that in the convention of Butcher and Cotter, the frequency arguments to 
the right of the semicolon may be writen in arbitrary order, hence we may in an equal description 
instead use 

or 

(3) f ^ 

for this description of the third order interaction between light and matter. 
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The polarization density 

The degeneracy factor K{—uj; u, ui, -co) is calculated as 

K{-io;uj,LJ, -u) = 2'+™-"p = 2i+°-33 = 3/4. 

From the reduced set of nonzero susceptibilities for the beam polarized in the xy-plane, and by using 

the calculated value of the degeneracy factor in the convention of Butcher and Cotter, we hence have 
the third order electric polarization density at Wo- = w given as p(")(r,t) = Re[PL"^exp(-iwt)], 
with 

p(3) = Ve (P^^h 

ft 

= {Using the convention of Butcher and Cotter} 

fj. L Q /3 7 

= {Evaluate the sums over (a;, y, z) for field polarized in the xy plane} 

Q 

= ,r„_/e fv^^) E^E'^E" +v(^) E^E^E^ + v^^) E^E^^E^ + v^^) E^^E^E^ 1 

^ ^'^^\^xxxx^bj^uj^-oj ^ Xxyyx-^u^u^-ui ^ Axyxy^u-^w^-w ^ A-xxyy^ui-^u^-ul 

+ e fv^^^ E^E^E^ +v(^' E^E'^E^ + v^^' E^E'^E'^ + v^^^ E^E^^E^ ]\ 

^ ^ylXyyyy-^w-^ul-^-w ^ A.yxxy^bJ^UJ^-UJ ^ A.yxyx^OJ^UJ^ -Ul ^ Xyyxx-^u-^u^-bJlf 

= {Make use of = E* and relations Xxxyy = Xy^L' etc.} 

= ,r„_/efv(3) R^IR^P -L v(3) £12/2™* , (3) |C;3/|2 , ^(3) rxiz-ij/pl 
^\^x\.Axxxx-^oj\-^u\ ^ Xxyyx-^U! ^ui ' A.xyxy\^ijj\ -^oj ^ Xxxyy-^oj\-^w\ J 

+ e fv^^^ E^li?"!^ + v'^^ 7-.J/* , (3) irsp™, (3) e;!/| e;^:^!-. 

~ ^yl/(.xxxx-^u>\-^uj\ ~ A.xyyx^w ^u> ~ A.xyxy\-^w\ -^ui ~ /ixxyy^u>\-^u>\ IS 

= {Make use of the intrinsic permutation symmetry} 

o 

^\.^xi\Axxxx\ uj\ ' Axxyyl oj\ / cj ' \Axxxx Axxyy J oj oj 
^y[iXxxxx\E^\ '^Xxxyyl^ujl ) E^ + {Xxxxx ^ "^Xxxyy) E^ . 

For the optical field being linearly polarized, say in the x-direction, the expression for the polar- 
ization density is significantly simplified, to yield 

P^f^So{3/^)e^/^ljEZ\'EZ, 

i. e. taking a form that can be interpreted as an intensity-dependent (~ l-E^P) contribution to the 
refractive index (cf. Butcher and Cotter §6.3.1). 

The wave equation - Time independent case 

In this example, we consider continuous wave propagation^ in optical Kerr-media, using light 
polarized in the x-direction and propagating along the positive direction of the ^;-axis, 

E(r, t) = Re[Et^(^;) exp(-ia;t)], ^u{z) = A^{z) exp{ikz) = exA^{z) exp{ikz), 

where, as previously, k = ujuq/c. From material handed out during the third lecture (notes on the 
Butcher and Cotter convention), the nonlinear polarization density for a;-polarized light is given 
asPr)=Pif\with 

P(f) = £o(3/4)e,xEL(-c^;o.,'^, -^)\EZ\''EZ 
= £o(3/4)xEL|E^|2e^ 
= eo(3/4)xEL|Ac.pA„exp(iA;2;), 



^ That is to say. a time independent problem with the temporal envelope of the electrical field 
being constant in time. 
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and the time independent wave equation for the field envelope A^^,, using Eq. (6), becomes 

^A(^ = ^^-^ fo(3/4)xEL|A^PA^exp(ifcz) exp(-ifcz) 



(3) I A |2a 

8c^fc ^^^^ 
■ {since k = u)no{cij)/c} 

i—yi^) lA Pa 



or, equivalently, in its scalar form 



Boundary conditions - Time independent case 

For this special case of unidirectional wave propagation, the boundary condition is simply a known 
optical field at the input, 

A^(0) = {known}. 

Solving the wave equation - Time independent case 

If the medium of interest now is analyzed at an angular frequency far from any resonance, we may 
look for solutions to this equation with |A„(2;)| being constant (for a lossless medium). For such 
a case it is straightforward to integrate the final wave equation to yield the general solution 

A^{z) = A„(0)exp[z^xilJA„(0)|2z], 

or, again equivalently, in the scalar form 

which hence gives the solution for the real-valued electric field E(r, t) as 

E(r,t) = Re[F,^{z) exp(-iwt)] 

= Re{At^(2;) exp[i{kz — u)t)]} 

= Re{AUO) eMii'^z + ^xilJ^S(0)|'^ -cot)]}. 

^ V ' 

= fectfZ 

From this solution, one immediately finds that the wave propagates with an effective propagation 
constant 

fceff = ^[no + ^xELK(0)|'], 
c ono 

that is to say, experiencing the intensity dependent refractive index 

"eflf = no + ^xELK(0)|' 
= no + n2|^S(0)p, 

with 
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